Spin-filtered edge states in graphene 
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Spin orbit coupling changes graphene, in principle, into a two-dimensional topological insulator, 
also known as quantum spin Hall insulator. One of the expected consequences is the existence of 
spin-filtered edge states that carry dissipationless spin currents and undergo no back-scattering in 
the presence of non-magnetic disorder, leading to quantization of conductance. Whereas, due to 
the small size of spin orbit coupling in graphene, the experimental observation of these remarkable 
predictions is unlikely, the theoretical understanding of these spin-filtered states is shedding light 
on the electronic properties of edge states in other two-dimensional quantum spin Hall insulators. 
Here we review the effect of a variety of perturbations, like curvature, disorder, edge reconstruction, 
edge crystallographic orientation, and Coulomb interactions on the electronic properties of these 
spin filtered states. 
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I. INTRODUCTION 

Since the fabrication of high mobility graphene based 
transistors, back in 2004^, the exploration of graphene 
has led to a number of striking experimental finding^ 
such as the quantum Hall effect^^^, even at room 
temperature, the Klein tunneling, the universal optical 
absorption, and also to an even larger amount of the- 
ory work predicting exotic electronic phases in graphene. 
Whereas most of the new electronic phases in Condensed 
Matter Physics are driven by Coulomb interactions, in 
2005 Kane and Mele (KM) predicted that spin orbit cou- 
pling (SOC) alone would turn graphene into a completely 
new electronic state, the so-called quantum spin Hall in- 
sulator (QSHI)^^^. Inspection of the energy bands of 
graphene, as described within the Kane-Mele model, look 
deceptively like an ordinary narrow gap semiconductor. 
However, a careful study of the wave function topology in 
reciprocal space reveals that the gap opened at the Dirac 
points by the SOC has different signs at the two val- 
leys. The twisting in the reciprocal state makes graphene 
topologically different from vacuum or graphene without 
SOC. The most conspicuous consequence of the topo- 
logically non-trivial band structure of the QSHI is the 
existence of spin-filtered edge states at the boundaries of 
graphene with other topologically trivial systems. 

The quantum spin Hall phase proposed by Kane and 
Mele can be portrayed as two copies of a quantum Hall 
phase^, one for each spin orientation and with opposite 
spin-dependent field orientations. Like in the conven- 
tional quantum Hall phase, the bulk states are gapped, 
but the edges hold states at the Fermi energy that carry 
current. In the quantum spin Hall phase, at a given edge, 
spin up electrons flow in the opposite direction to that of 



spin down electrons, resulting in a pure spin current flow. 
Since backscattering in a given edge requires reversal of 
the spin, the conductance of the edge states is expected 
to be quantized even in the presence of time-reversal pre- 
serving disorder. 

In their original work, KM studied the competition be- 
tween the intrinsic SOC, which drives the QSHI phase, 
and inversion symmetry breaking perturbations, such as 
a Rashba spin orbit term and a staggered perturbation, 
which make graphene an ordinary insulator. The in- 
terplay of these couplings with non-magnetic disorder 
was later addressed by L.Sheng et al.^ and D.N.Sheng et 
al.^ . They found that, as expected, the QSHI phase, and 
thereby the spin- filtered edge states, were robust with 
respect to the addition of time-reversal symmetric disor- 
der. 

Here we address four different topics: first, the rela- 
tion of the KM model with the microscopic description 
of SOC in graphene. In section II, we show that the KM 
Hamiltonian, originally derived from heuristic symmetry 
considerations, is a good effective one orbital description 
of the low energy bands of flat graphene, as described 
with a 4-orbital Slater-Koster (SK) Hamiltonian includ- 
ing SOG^-i^ 

Our second topic refers to the influence of the geometry 
of the sample on the the electronic properties of the spin- 
filtered edge states. We consider the effect of the crys- 
tallographic orientation, either zigzag or armchair, the 
effect of a reconstruction into the so-called reczag edge, 
and the effect of the curvature. In the original papers, 
KM studied the spin-filtered states in graphene termi- 
nated by zigzag boundaries. As we discuss in sections Hill 
and[Vl the choice of zigzag ribbons with pre-existing non- 
spin filtered edge states, obscures somewhat the dramatic 



effect of SOC on the electronic properties of graphene. 

The third topic addressed in section llVl is related to the 
robustness of the edge state conductance with respect to 
the presence of non-magnetic disorder. We study the ef- 
fect of various types of disorder, Anderson type and edge 
roughness, on the two terminal conductance of zigzag 
terminated ribbons. 

Finally, our fourth topic concerns the study of the ro- 
bustness of the edge states with respect to magnetic or- 
der. In the case of zigzag ribbons ignoring SOC, the 
edge states undergo a Stoner inestability^^^^ for arbi- 
trary small electronic repulsion. Thus, the robustness 
of the spin- filtered edge states is compromised^^, unlike 
in the case of armchair terminations. In section IV we 
discuss the interplay between SOC and electron-electron 
interactions in both geometries. 



II. SPIN ORBIT IN GRAPHENE AND THE 
KANE MELE HAMILTONIAN 

The KM Hamiltonian can be seen from two different 
perspectives: as a theoretical model that describes a new 
phase of matter, the QSHI phase, or as a effective model 
Hamiltonian to describe the effect of SOC on the tt elec- 
trons of graphene. The KM model is inspired by a mathe- 
matical model for spinless fermions under the influence of 
an inhomogeneous magnetic field that, with a vanishing 
total flux, gave rise to a quantized Hall response. This 
model, originally proposed by Haldane^^, breaks time re- 
versal invariance but shows a way to obtain quantized 
Hall effect without an applied external field. By including 
the spin degree of freedom and giving to the inhomoge- 
neous magnetic field a spin dependent value, KM found 
a way to restore time-reversal symmetry in the model, 
which would yield the so called QSHI phase. 



A. Spin orbit coupling in graphene 

We briefly discuss the fundamentals of SOC in 
graphene. The low energy physics of graphene, ignoring 
SOC, is described by the tt orbitals, with orbital angular 
momentum / = 1 and m = 0, if we chose the quantization 
axis perpendicular to the plane. The SOC, as described 
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FIG. 1: (Color online) Left panel: Graphene band structure 
computed with the 4-orbital orthogonal tight-binding model 
with a SOC of A=4eV (blue lines) and without SOC (green 
lines). Dashed and solid lines represent empty and occupied 
bands, respectively. Right panel: zoom of the band structure 
at K point. 



carbon atom and a total of 16 in the unit cell of the 
honeycomb lattice. In this case, the 16 x 16 Hamilto- 
nian can be obtained using a SK parametrization of the 
tight-binding Hamiltonian^iiSii^. In figure 1 (left panel) 
we show the energy bands of graphene so obtained. For 
the sake of clarity we use an unrealistically large SOC, 
A = 4eV, 3 orders of magnitude above the accepted value 
A c^ lOmeV. Whereas the large energy scale is not much 
affected by the SOC, a band gap opens right at the K 
and K' points (right panel) . For realistic values of the 
SOC, A c^ lOmeV^ the value of the gap is in the range of 
lO/j^eV. Thus, from a practical perspective, the SOC gap 
in graphene can be ignored in most instances. 

Within this 4-orbital picture, the smallness of the band 
gap comes from the fact that A is small and it only affects 
the TT bands to second order^'^^. Although it does not 
make a difference in the rest of the discussion, it is worth 
pointing out that the 4-orbital picture is not complete 
because the d orbitals happen to play a role. The xz 
and yz orbitals are weakly hybridized to the tt orbitals, 
even for zero SOC. The d orbitals, in turn, have their 
own SOC, as a result of which the ir-d mixing is spin 
dependent. This results in a contribution to the tt bands 
with the same symmetry than the one coming from the A^ 
couphng through the a bands, but twenty times higher—. 



couples, to lowest order in A, the tt electrons with spin 
s^ to / = 1,771 = ±1 electrons with opposite spin Sz- 
Thus, SOC mixes the tt electron with the a electrons. A 
possible approach is to give up an effective description 
in terms of the tt orbitals only and use a tight-binding 
model that includes the four atomic orbitals, i.e., the 
three responsible for the sp"^ hybridization and the tt or- 
bital, in each atom. Thus, there are 8 spin orbitals per 



B. The Kane Mele Hamiltonian 

The results above suggest that some sort of effective 
Hamiltonian where the higher energy bands are inte- 
grated out can describe the low energy physics around 
the K points in terms of the tt orbitals only. The numer- 
ical calculations and symmetry considerations show that 



close to the K points Sz is a good quantum number. In 
the lattice representation, the effective model that does 
the job was proposed by KM: 



Reczac zz(57) 



HkM = tKM /Z ^i^^^f^J ^ ^ik)Cj 

{{hJ)) 



(2) 



where the SOC yields a spin-dependent term which con- 
nects second neighbors in the honeycomb lattice through 
the two vectors dkj and dik as shown in o figure ([2j), and 
txM is the spin-orbit coupling strength. When summed 
to the standard first neighbor hopping t term, the mo- 
mentum representation of the Hamiltonian for graphene, 
including SOC, can be written as: 



n{k) = tKMg{k)sz(Tz + tf{k)ax 



(3) 



where ax and a^ are the Pauli matrices in the sub lattice 
space, t is the first neighbor hopping, f{k) = 1 -\- Exp[ik- 
ai] + Exp[ik ■ a2], g{k) = i{—Exp[ik ■ ai] + Exp[ik - a2] — 
Exp[ik • as] + Exp[ik • ai] — Exp[ik • a2] + Exp[ik • as]) 
and txM is the spin-orbit coupling strength. 

Notice that g{k) = —g{—k) because the sign of the 
spin-dependent second-neighbor hopping is also sensi- 
tive to the clockwise-anticlockwise nature of the hopping 
path. At the Dirac points f{k) vanishes and the effect 
of the SOC Hamiltonian is to open a valley and spin de- 
pendent gap given by the following Hamiltonian: 



HsO = ^SO'^zO'zSz 



(4) 



where A so = Sa/S^km is the band gap divided by 2, r^ = 
±1 for i^ and K^ valleys, a^ = ±1 for A and B subltatice 
and Sz = ±1/2 for the spin projection perpendicular to 
the plane of graphene. At the K and K^ valley the top of 
the valence band and bottom of the conduction band are 
sub lattice polarized. It is quite apparent that, for a given 
spin orientation Sz = +^, the wave functions change sub 
lattice character as we change the valley index. Thus, 
the gap changes sign, at a given spin, when we go from 
K to K' . This twist in the reciprocal space makes the 
SOC gap topologically different from the gap opened by 
a sublattice symmetry breaking or staggered potential 



n 



stg 



Atz which favors the same sublattice at both 



valleys. 



III. EDGE STATES: THE EFFECT OF 
CRYSTALLOGRAPHIC ORIENTATION 

Whereas the QSHI phase is a bulk property, it is re- 
vealed at the boundaries through the appearance of edge 
states, very much like it happens in quantum Hall bars^^. 
In most instances edges states are studied by looking at 
sufficiently wide strips with two edges rather than semi- 
infinite systems with only one edge. Strips or ribbons 
are one-dimensional crystals whose electronic structure 
is simplified taking advantage of the translational invari- 
ance along the longitudinal (edge) direction, which we 




Zigzag 

FIG. 2: (Color online). The three types of edges consid- 
ered: reczag (top), armchair (sides), and zigzag (bottom). 
The bond-vectors di featured in the Kane-Mele Hamiltonian 
are also shown. 



take to be along the x axis. The existence of edge states 
is an expected consequence of the non trivial topological 
order of the QSHI phase. However, the study of their 
properties in ribbons of the honeycomb lattice has some 
subtleties associated both with inter-edge coupling, on 
account of their finite transverse size, and with the ex- 
istence, along some crystallographic directions, of edge 
states even in the absence of SOC^^. Here we discuss 
the edge states in three different types of termination: 
zigzag, armchair and reczag, using the two different tight- 
binding models described above, namely, the effective 1- 
orbital KM and the 4-orbital SK model. We also study 
the effect of curvature on the electronic properties of the 
zigzag edges. 



A. Zigzag edge states 

The most studied graphene termination is the so-called 
zigzag edge [shown in Fig. (|2])], which presents pre- 
existing edge states, i.e., edge states even without SOC. 
This can be seen in Fig. [3l where we plot the band struc- 
ture of a zigzag ribbon with A^ = 48 sites in each unit cell 
calculated within the KM model (identical results are ob- 
tained within the SK model"*^^). We consider two cases, 
one without SOC and one with a finite SOC . In both 
cases the bands have a two-fold degeneracy associated 
with the spin. (In the case without SOC the degeneracy 
associated with the existence of two identical edges is 
slightly removed by the interaction between edges.) The 
wave functions (modulus square) across the y direction 
of the highest energy valence band as a function of the 
longitudinal momentum k are shown in the lower pan- 
els. Without SOC, the flat mid gap bands correspond 
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FIG. 3: (Color online). Band structure of a graphene zigzag 
ribbon for (a) txM /t — and (b) txM /t — 0.03. The unit 
cell contains N — 4& atoms. Panels (c) and (d) represent 
the highest valence band charge density |^]^^|^ for each site i 
across the ribbon width. Red and blue colors refer to different 
spin projections. When the SOC is turned on, the valence and 
conduction bands acquire some dispersion between K and K' 
points and become spin- filtered. 



to localized states at the edges^^. In the case of zigzag 
ribbons, the valley character is preserved^. Thus, it can 
be said that edge states go from K to K' in momentum 
space. In the presence of SOC (right panels) the valence 
and conduction bands connecting K and K' points in 
the Brillouin zone acquire an apparently innocent dis- 
persion, but the wave functions undergoes a dramatic 
change. Whereas in both cases the states at the valence 
and conduction bands along the path K ^ K' are lo- 
calized at the zigzag edges, in the presence of SOC they 
become spin- filtered, giving rise to helical edge states^^^. 
This means that electrons with opposite (same) spin carry 
current in opposite direction along the same (opposite) 
edge. Interestingly, the occupied bulk states also carry 
a spin current in the opposite direction so that the net 
spin current at the edge is zero^^ 




FIG. 4: {Color online) Band structure and squared highest 
valence band wavefunction |^]^^|^ for an armchair graphene 
ribbon with A/" = 50 atoms across the ribbon width. Red and 
blue colors refer to different spin projections. For panels (a) 

and (c), txM /t = and for (b) and (d), txM /t = 0.1. 



width^^^^^. In the lower panels we show the wave func- 
tions (modulus square) for each spin projection of the 
highest valence band with no SOC. It is apparent that 
they are delocalized along the transverse direction. 

By making the ribbon sufficiently wide, the confine- 
ment gap is smaller than the SOC gap. When this hap- 
pens, two bands crossing the Fermi energy appear whose 
wave functions are localized at the edges. In addition, 
their wave functions [Fig. lU^d)] show the same correlation 
between velocity, spin orientation and edge than their 
counterparts in the zigzag ribbons. Thus, spin-filtered 
edge states appear without having preformed edge states. 
Similar results have been obtained by E. Prada et al?^^ 
using the kp approach, originally devised for ribbons 
without SOC^^, complemented with the SOC term in 
Eq. dH). 



C. Reczag edge states 



B. Armchair edge states 

We now consider a different crystallographic orienta- 
tion of a honeycomb lattice, the armchair edges [see Fig. 
[2]. In Fig. Hlwe plot the energy bands of the KM model 
for a semiconducting armchair ribbon with A/" = 50 sites 
in the unit cell, both with null (left) and finite (right) 
SOC. In contrast to the zigzag ribbons, this termina- 
tion lacks pre-existing edge states. Within this model 
and with zero SOC, the band structure of armchair rib- 
bons is either conducting or insulating depending on the 



We now consider reconstructed edge states in zigzag 
graphene ribbons. This is motivated by density func- 
tional theory calculations claiming that zigzag edges are 
prone to edge reconstructions^^. Here we consider the 
so-called reczag, which is a self-passivated zigzag edge, 
where the edge hexagons change into a configuration of 
alternating pentagon and heptagons. In order to study 
the electronic structure and the edge states of these sys- 
tems, we make use of our 4-orbital orthogonal SK model 
with SOC, taking into account that the hopping depends 
on the distance^^. This is the a necessary methodology 
here because a the 1-orbital KM model cannot account 



for the physics of the newly formed atomic structure. 
Density functional theory calculations including SOC are 
also valid, but may require in some cases great computa- 
tional resources. 
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FIG. 5: {Color online) Band structure and squared highest 
valence band edge states wavefunctions |^J^^|^ for a reczac 
graphene ribbon with 104 carbon atoms per unit cell. For 
panels (a) and (c), A =OeV and for (b) and (d), A = 4eV. 



We obtain the geometrical disposition of the atoms 
from Ref. [23 and the SK parameters from Ref. [l9|. 
The electronic configuration for the low energy spectra 
is shown in the Fig. [5] both with null (left) and finite 
(right) SOC. Because of the reconstruction, the unit cell 
of the reczag ribbons is twice as large along the longi- 
tudinal direction and, thereby, the Brillouin zone is half 
as large as the one for the zigzag ribbons. As a result, 
the edge states occupy a larger percentage of the Bril- 
louin zone. It is apparent that the band edges are no 
longer electron- hole symmetric. This is due to the loss of 
the bipartite character of the reconstructed lattice. This 
also results in dispersive edge states, in contrast to the 
zigzag case. In the lower panels of Fig. [5] we show the 
wave function map for the lowest energy conduction band 
states projected on each atom as function of k. Without 
SOC, the reczag ribbon presents edge states that are not 
spin-filtered and present a four-fold degeneracy (edge and 
spin). When the SOC is turned on, the four- fold degen- 
eracy is lifted and they acquire a spin-filtered character, 
shown in Fig. [5jd). In this instance, the non-monotonic 
character of the band velocity makes the effect more in- 
tricate. As a consequence of the shape of the bands, each 
of the two spin-filtered edge states cuts three times the 
Fermi level, introducing more conduction channels near 
the Fermi energy. 



D. The effect of curvature 

The effect of SOC in graphene is small due, in part, to 
the decoupling of the tt and a bands for fiat graphene. 
Since it is possible to have curved graphene ribbons, 
made, e.g., by opening carbon nanotubes^^^^^, the in- 
terest on the effect of curvature on the spin filtered edge 
states has increased^^. In a curved ribbon, the tt orbitals 
are not longer decoupled from the a orbitals. This is 
similar to the case of fiat graphene under the action of a 
perpendicular electric field^. As a consequence, we can 
expect a new terms in the effective hamiltonian for the tt 
linear in the atomic SOC constant A, as a result of which 
Sz is not a good quantum number anymore. In the case 
of fiat graphene under an external electric field the new 
term is the so-called extrinsic SOC or Rashba interaction, 
which acts as a spin-fiipping nearest-neighbor hopping. 

The electronic structure of curved graphene ribbons 
with zigzag terminations is shown in in Fig. [6] and 
has been obtain within the 4-orbital SK model includ- 
ing SOGi^. The low energy band structure (not shown) 
remains similar to the case of fiat ribbons, except for a 
the appearance of a small dispersion in the originally fiat 
bands. For finite SOC, the edge bands are spin filtered 
as well, but with an important difference with respect to 
the flat case. Since Sz is not a good quantum number, 
the spin polarization of these states acquire an angle, ^, 
with respect to the normal direction at each atom of the 
edge. In Fig. [6fa) we represent the wavefunction of the 
edge states at /c = ±0.997r/a. The fact that the spin 
polarization points in opposite direction for these two k 
states reveals they form a time reversal Kramers pairs. 
In Fig. [6](b) we also show how 6 evolves with respect 
to the curvature for two different values of the SOC. We 
observe that the effect of curvature is dramatic for a real- 
istic value of the spin-orbit coupling, where almost neg- 
ligible curvatures make the spin point in the tangential 
direction. 

As a conclusion for this section, we have seen how SOC 
gives rise to the appearance of spin filtered edge states, 
regardless of the crystallographic direction of the edge. 
Interestingly, the bands of the edge states do depend on 
the direction. Unfortunately, the energy scale for the 
edge states dispersion is below ImeV when computed 
with realistic values of the SOC for graphene. However, 
understanding spin filtered edge states in this chemically 
simple situation could shed some light on the edge states 
in more complicated systems. 



IV. DISORDER AND ELECTRONIC 
TRANSPORT 

The topological character of the QSHI phase manifests 
itself in the robustness of the spin-filtered edge states. In 
the presence of non-magnetic disorder they behave like 
a pair of robust quantum channels, presenting a conduc- 
tance quantization. Go = 2e^//i, even when Sz is not a 
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FIG. 6: (Color online) (a) Spin polarization of edge states 
at /c = 0.997r/a(top) and k = -0.997r/a (bottom), (b) Angle 
between the normal direction of the edge and the spin po- 
larization as a function of the curvature for A = 5meV (red 
dots) and A = 2eV (blue squares). 



good quantum number^^^. Shortly after KM proposed 
the existence of a QSHI phase in graphene, the conduc- 
tance quantization was numerically verified by L. Sheng 
et at. for the Anderson disorder model^*^. Here we carry 
out a similar study through an analysis of the transport 
properties of zigzag ribbons in a variety of situations. We 
are mostly intrigued by the case of curved ribbons as the 
ones introduced in previous section. First, for illustra- 
tion purposes, we revisit the problem in fiat ribbons. As 
shown in SecHIl zigzag graphene ribbons with SOC con- 
stitute the most simple model of QSHI since they present 
a single branch of helical states for any width of the rib- 
bon and only inter-edge backscattering could affect the 
quantization of conductance. Finally, we ask ourselves if 
the spin non-conserving effect of the curvature can play 
any role in the inter-edge backscattering. 

We study two types of disorder: Anderson type, where 
we sum random on-site values e^ distributed uniformly 
in a range [—W^ W] to the Kane-Mele Hamiltonian, and 
structural defects, such as constrictions across the ribbon 
width. We are interested in coherent transport so we 
make use of the standard Landauer formalism as in Ref. 
|3Q| : The disordered region contains N = Lx x Ly = 300 
lattice sites, where L^y = 10 is the number of horizontal 
lines along the ribbon width and L^ = 30 the number of 
atoms along each chain. 

In Figure [Tl^a) we show the density of states at the 
Fermi level in real space for increasing disorder strength 
W and for a SOC strength ^km = O.lt. All results with 
W y^ have been averaged over 10 random realizations 
of disorder. In the trivial insulator case for ^km = 0, 
where a bulk gap arises due to lateral confinement, dis- 
order induces intra-edge backscattering, localizing the 
edge states and strongly reducing the conductance at the 
Fermi level [see Fig|7fb)]. The effect of SOC on the con- 
ductance (close to the Fermi energy) is dramatic: even 
for large values of W^ backscattering is almost entirely 
suppressed, as expected from their spin-filtered charac- 
ter. 
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FIG. 7: {Color online) (a) Density of states at the Fermi en- 
ergy projected on real space for a zigzag graphene nanoribbon 
with tKM = O.lt and increasing values of the disorder strength 
(top) W = 0.0, (middle) W = 1.0 and (bottom) W = 2.0. (b) 
Conductance at the Fermi energy with (bottom) and without 
(top panel) SOC for increasing disorder strength W . Error 
bars in conductance reflect statistical dispersion over disorder 
realizations. 



In finite width systems, when spin-filtered states of dif- 
ferent edges having the same spin polarization are able to 
interact one with each other, a finite inter-edge backscat- 
tering will be present and an exact quantization is no 
longer expected. We study this effect by introducing con- 
strictions in the ribbon (see Fig. [51 right panels). We 
model them by removing atoms in an area near the edge 

2 

delimited by gaussian functions /3e~^^ {x runs along the 
edge) for random /3 and a. Here we take ^km = O.Olt 
or A = 4 eV. In Fig. [8] we show the conductance for 
three different constrictions as a function of energy. We 
can observe that, while in the absence of SOC (dashed- 
dotted lines) the conductance is negligible at the Fermi 
energy because the edge states are able to interact with 
the states in the opposite edge and backscatter, for finite 
SOC the conductance increases appreciably (solid lines). 
These results have been obtained within the KM model. 

L. Sheng et alJ numerically showed that the QSHI 
phase survive in presence of non conserving spin terms. 
Here we study this by introducing a transversal bending 
over the constrictions previously analyzed. The curva- 
ture introduces a non trivial angle for the spin polariza- 
tion of the spin-filtered edge states of graphene nanorib- 
bons as mentioned in section IIIII These angles are dif- 
ferent for each edge. This behavior raises the question 
whether the relative phase between the spin orientation 
of the two channels could affect the process of backscat- 
tering. In order to answer this question we need to use 
our 4-orbital SK Hamiltonian instead of the single or- 
bital KM model. First, we analyze the transport for flat 
constrictions and compare the results with the KM ones. 
To be meaningful we need to saturate the a dangling 
bonds of the carbon atoms with hydrogen atoms. We 



compute the two terminal conductance with the ANT ID 
code, which is part of the ALACANT package^^. In Fig. 
[Hl^a-c) the dotted Hues show the conductance for the three 
constrictions considered. The difference between the SK 
and the KM conductance reveals the approximation in- 
herent to the KM model, but are essentially identical. 
While the SOC term is exact in the SK formalism, the 
SOC introduced in the KM Hamiltonian is a second or- 
der perturbative approach to the exact value which still 
preserves electron- hole symmetry. ^^. 

Finally we consider the same constrictions with a 
transverse bending, with a radii of curvature of 4.1nm 
(see right panels in Fig. [8]). In this situation we ex- 
pect that the backscattering is less effective due to the 
fact that the spin polarization are not parallel between 
the channels of opposite edges. The results are shown 
by dotted and dashed lines, revealing that the process of 
backscattering is the same for both situations, namely, 
for a flat ribbon where Sz is a good quantum number and 
for a curved ribbon where Rashba like terms are present. 
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FIG. 8: {Color online) (a-c) Conductance for three different 
types of constrictions computed in different situations. Solid 
and dashed-dotted lines correspond to flat constrictions using 
the KM Hamiltonian with and without SOC, respectively. 
Dotted lines show the conductance for the same constrictions 
computed using SK approximation. Dashed lines are related 
to the same cases but with transverse bending as shown in 
panels (d-f), using the SK approximation. 



V. THE EFFECT OF COULOMB 
INTERACTIONS 

In the previous section we have seen how edge states 
in the Kane-Mele model are robust with respect to non- 
magnetic weak disorder. In this section we address 
if, as predicted^^^, they are also also robust with re- 
spect to Coulomb interaction. For that matter, we 
use the Hubbard model within a mean field approxima- 
tion (MFA), where the interaction term is approximated 
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FIG. 9: {Color online) On the left panel we show the 
phase diagram of an AC ribbon with Ny = 50 in presence 
of Coulomb {U/t) and spin-orbit {txM /t) interactions. The 
labeled regions are related to the following phases: (I) Insulat- 
ing or semiconducting phase, (II) antiferromagnetic insulating 
phase, (HI) Quantum spin Hall insulating phase with gapless 
edge states and (IV) antiferromagnetic phase where metallic 
edge states are still preserved. On the right panel we show 
the phase diagram for the AF insulating (AF), the AF valley 
half-metal (HM) and the non-magnetic (NM) phases for the 
ZZ ribbon with Ny = 24. 



by U^^{ni^^{ni^) -^ ni^^{ni^^). More sophisticated ap- 
proaches have also implemented recently^. 

At half-filling, the Hubbard model in the honeycomb 
lattice has magnetic order at sufficiently large U . In the 
MFA Uc = 2.2t, using Quantum Monte Carlo Uc > 5t^^. 
This magnetic ordering is always accompanied by a gap 
opening transforming graphene into a Mott insulator. On 
the contrary, the zigzag edge states have a very large 
density of states at the Fermi energy so that, at zero 
SOC, undergo the ferromagnetic transition at arbitrary 
small U. In the case of finite width ribbons the inter- 
edge correlations are anti-ferromagnetic and a band-gap 
opens in the otherwise conducting band structure^^^^. 

Thus, under the separate influence of SOC and 
Coulomb repulsion, the preformed zigzag edge states un- 
dergo a radical change in their electronic structure. In 
the first case they carry spin-currents and acquire a linear 
dispersion, in the second case they are ferromagnetic and, 
for finite width ribbons, open a gap. Thus, it was not 
clear a priori what would be the combined effect of SOC 
and Coulomb repulsions. This problem was addressed by 
two of us^^. When the SOC is turned on, the effect on 
the insulting band structure of the magnetic phase is to 
decrease the band-gap in one valley and increase it in the 
other. This valley symmetry breaking goes together with 
a reduction of the edge magnetic moment. As the SOC 
is increased, the gap closes in on of the valleys, the edges 
being still magnetic. This phase is a valley half-metal. 
As the SOC increases further, the magnetic moment goes 
to zero and the non-magnetic phase is recovered. These 
findings are summarized in the phase diagram of figure 
(|9j)b, for a ZZ graphene ribbon containing Ny = 24, that 
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FIG. 10: (Color online) Band structure and projected spin 
density along the valence band of an AC ribbon with Ny — 
50 for three different states along the line Ikm = O.lt in 
the phase diagram (a) (U = l.Ot,tKM = O.lt), (b) (U = 
3.0t,tKM = O.lt) and (c) (U = 3.5t,tKM = O.lt) 



shows that, because of the SOC, a finite U is needed to 
make the edges magnetic. 

We now consider the armchair termination, for which 
no preformed edges states exist, within the MFA for 
the Hubbard-Kane-Mele model. This system has also 
been addressed using bosonization techniques^. In fig- 
ure ([T0|) we show our results for the the band structure 
and edge-state wave function of an AC terminated rib- 
bon, for 3 values of /7, and txM = O.lt. At the left, with 
U = t^ the system is non-magnetic, and very similar to 
the U = case. When the Coulomb interaction is in- 
creased above a critical value, the bulk atoms undergo 
an antiferromagnetic transition but the edge atoms stay 
non-magnetic and the spin-filtered edge states are pre- 
served. Only when the interaction is increased above a 
Uc significantly larger than the one that makes the bulk 
magnetic, the edge atoms go magnetic as well and the 
spin- filtered edge states dissappear. Thus, for AC edges 
it is fair to stay that spin-filter edge states are robust 



with respect to Coulomb repulsion, since they survive 
even when the bulk goes antiferromagnetic. 

In figure ([9|)a we show the phase diagram of a semicon- 
ducting AC ribbon in presence of both Coulomb and spin- 
orbit interaction. For txM < 0.05t and U < Uc (region I) 
the ribbon remains semiconducting and for U > Uc the 
gap increases and the ribbon becomes an antiferromag- 
net (region II). For Ikm > 0.05t the ribbon transform 
into a QSHI holding spin filtered edge states as those 
shown in the case of the ZZ ribbons (region III) and, as 
the Coulomb interaction increases, we find a new phase 
(region IV) which is antiferromagnetic with a magnetic 
moment per atom in the edge of \m\ ~ 0.2 but where the 
gap is still closed and the states remains spin- filtered. 
VI. SUMMARY 

We have studied the electronic properties of the spin- 
filtered edge states that appear in graphene because of 
spin-orbit coupling. Whereas they are a generic feature 
of the Quantum Spin Hall Insulator phase, their prop- 
erties are very different depending on the existence of 
preformed edge states, i.e., edge states that occur even 
with no SOC. For instance, we reczag edges have 3 con- 
ductance channels, and zigzag edges are not robust with 
respect to Coulomb interactions. Whereas graphene is 
not a good material to study spin-filtered edge states ex- 
perimentally, because of the low spin orbit coupling, it 
provides an ideal platform to study generic features of 
such an interesting object which can be found in other 
two dimensional topological insulators. 
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